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DESCRIPTION 

METHOD OF POSITION DETERMINATION 

The present invention relates to a method of position determination 
in a radio system and to apparatus for position determination in a radio 
5 system. 

European Patent Publication 0 654 677 A describes a method of 
processing a signal in a multi-path environment. 

However, such methods may not perform particularly well In the 
10 presence of diffuse reflections especially when such components are very 
plentiful and closely spaced, as for example occur in indoor environments 
involving rough uneven surfaces and clutter. Accordingly, imprecise 
estimations and hence inaccurate positioning result. 

An object of the present invention is the avoidance or elimination of 
15 diffuse multi-path inaccuracies. 

Another object of the present invention is to provide good 
positioning accuracy. 

Another object of the present invention is to remove diffuse 
scattering components from received signals. 

20 

According to the present invention, there is provided a method of 
position detemiination in a radio system, the method comprising 
correlating a signal received at a unit with a replica signal at the unit, and 
processing the correlated signal with an optimisation function comprising 
25 an exponential term in combination with a second term. 

In this way, the method of the present invention may provide 
accurate position measurement in indoor environments. 

Preferably, the exponential term is in the form Be""^, and the second 
term is of the form: 
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The method may comprise effecting an integration with the repiica 
signal, and/or fitting the optimisation function and a Line-of Sight 
5 correlation function with a set of parameters, and/or superposing the 
diffuse correlation output with a Line-of-Sight function output and fitting 
with correlation data of laiown values for the Line-of-Sight output. 

The method may comprise first operating a multipath mitigation 
technique to effect correlation of the received and replica signals. 
10 The multipath mitigation technique may comprise a Multipath 

Estimating Delay Locks Loop (MEDLL) technique, or a Minimum Mean 
Square Enror (MMSE) technique. 

According to the present invention, there is also provided a 
computer program product direcHy. loadable into the Intemal memory of a 
15 digital computer, comprising software code portions for performing the 
method of the present invention when said product is run on a computer. 

According to the present invention, there is also provided a 
computer program directly loadable into the intemal memory of a digital 
computer, comprising software code portions for performing the method of 
20 the present invention when said program is run on a computer. 

According to the present invention, there is also provided a canier, 
which may comprise electronic signals, for a computer program 
embodying the present invention, and/or electronic distribution of a 
computer program product, or a computer program, or a carrier of the 
25 present invention. 

According to the present invention, there is also provided apparatus 
for position determination of a radio system, the apparatus comprising 
means to correlate a signal received at a unit with a replica signal at the 
unit, and means to process the correlated signal with an optimisation 
30 function comprising an exponential term in combination with a second 
term. 
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The present invention is applicable to radio systems especially 
involving GPS and GSM networks, particularly in indoor environments and 
can be Implemented in hardware and/or software. 

In order that the present invention may more readily be understood, 
5 a description is now given, by way of example only, reference being made 
to the accompanying drawings in which:- 

Figure 1 is a schematic diagram of a receiver unit of the present 
invention; 

Figure 2 is a hardware implementation of the unit of Figure 1 ; 
10 Figure 3 is a flow diagram of the signal processing operation of the 

method of the present invention; and 

Figure 4 depicts the surfaces that drffusers lie on. 

Figure 1 shows a mobile radio receiver unit 1 of the present 
15 invention which is located in a room 2 with a variety of objects including 
chairs 3, desks 4, bookcase 5 and reflective surfoces e.g. walls 6. These 
objects act in combination to generate a large number of diffuse reflections 
for GPS signals of the type for reception at unit 1 . . 

Details of the components in a hardware implementation of unit 1 
20 are shown in Figure 2, having an antenna 10 to receive incoming signal 
R(t) and pass it to multiplier block 11 for multiplication with output C from 
the carrier generator 12 in order to mix down signal R(t). The resultant 
signal is passed to a block of correlators 13 for sequential correlation of 
the mixed down signal and the replica code from code generator, in order 
25 to produce correlation data outputs from the I (imaginary) and Q (real) 
channels using input signals from code generator 14. 

The correlation data is then passed to signal processing block 15 
which adjusts appropriately both the code generator 13 and carrier 
generator 12. 

30 Figure 3 is a flow diagram of the major steps in the correlation 

operation for signal processing to remove the diffuse components as in the 
present invention. 
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Step 1 of the flow diagram involves standard signal processing for 
position measurement in the presence of multi-path signals, typically 
accomplished by the Muitipath Estimating Delay Lock Loop (i\^EDLL) 
method generally as described in European Patent Publication No. EP 0 
5 654 677 A. 

Step 1 involves using this method to provide an optimal fit using N 
triangles" thereby yielding estimates for the muitipath parameters:- 



Then step 2 adds to the standard correlation operation by using the 
initial value for t_0 obtained from the MEDLL method, incorporating small 
initial values for B and a, and then correlating the function 



15 with the replica signal s(t-r) . 

The correlation is performed using multiplier block 11 so that the 
received signal r(t) is multiplied by the replica code on the carrier wave to 
mix down the received signal as well as obtain the correlation data. This is 
then summed by integrator 16 in Step 3 to give correlation function data for 
20 the I and Q channels, the correlation function being complex due to the 
phase of the signal. 



10 



{a_0,....a_N} 

{t^1.....t^N} 

{q_1,....t^N} 




(1) 



In step 4, f (t) is fit with the Line-of-Sight triangle with parameters 
to the correlation data using known values 



25 



{a_1, . .a_N} 
{t^1.....t^N} 
{CU1....CLN} 



where 



30 



ao,...aN are the amplitudes 
To,....T;^r are the delays 
qo,....qN are the phases 
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of all the reflected components In the received signals. 
Step 5 involves finding values for 

B, 
a. 

5 a_0, 
clO. 

and E (the mean square error), 
trying different values until a global minimum error is reached, and 
10 storing all the derived values. 

Step 6 Involves Integrating the same operations for the respective 
values, and then feeding them back to Step 2 to reiterate the process. 

Thus the function f(T) which is characterised in equation (1) acts 
efficiently to clean the correlation functions of diffuse reflections. 
15 As a variant, a software implementation of the correlation operation 

incorporates an algorithm to effect: 

• performing the correlation of the local replica code with the received 
signal data using multiplication and summation of the sampled values 
of the local code and the received signal, to produce demodulated 

20 correlation data; 

• operating the conventional position measurement method by fitting N 
amount of triangles of fixed width 2* chip width to the demodulated 
correlation data such that N triangles give the minimum global error; 
and 

25 • operating the f (t) of equation 1 and Stage 2. 

In greater detail, Figure 3 represents an additional stage of signal 
processing, in accordance with the present invention, for the purpose of 
cleaning a received correlation function of remaining diffuse reflection 
components for a number of multipath mitigation techniques, for example a 

30 Multipath Estimating Delay Lock Loop (MEDLL) or a Minimum Mean 
Square Error (MMSE) technique. 
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Although the method shown In Rgure 3 and the following 
description both refer to a l\AEDLL technique, it should be understood that 
the improvement to this technique described hereinafter is equally 
applicable to other forms of multipath mitigation technique, one example of 
6 which is a technique referred to as the Minimum Mean Square Error 
(MMSE) mettiod. 

The additional stage in the signal processing utilises a cleaning 
function, /(r), in the maximum likelihood estimator (MEDLL) or Bayesian 
estimator (MMSE), making possible the removal of the diffuse scattered 
10 components from the received correlation function data. 

The present invention allows for the optimisation of the 
implementation of a maximum likelihood estimator, using a novel 
correlation function that is a sum of a line-of-sight (LOS) correlation 
function and a diffuse correlation function, both of which are dependent on 
15 the LOS delay Tq . 

The improved estimator utilises an approximation of the diffuse 
scattering process. 

The maximum likelihood estimator of a conventional MEDLL 
technique as currently in use is given by: 



20 




(3) 



where 




and where P(r„) is a triangle of unit height centred at time r =t„ 
25 with a base of two chip widths 2T^ , 



phase of the nth reflection 9^ 



nth amplitude a„ and 
nth delay r„ 
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Which are parameters in the correlation function modeL 
The term Z)(r), where: 

D(T)^jr(t)s(t^r)xIt . (5) 
results from the received signal data, r(0. correlated with the 
5 replica signal, s(t) , in the receiver to obtain the correlation function data. 

Thus, the maximum likelihood estimator for the multipath 
parameters {^o- -^jv^^o- -^/^^^o - •^/^} "sed in the MEDLL method Is given 
by those parameter values that minimise the estimator L . 

In deriving the f(r) function of equation (1), three main factors were 
10 considered, namely: 

i) the combined amplitude of randomly-scattered superposed rays as 
path length or time increases; 

ii) the number of randomly-placed diffusers as path length or time 
increases; and 

15 iii) the free-space decay of ray amplitudes as path length or time 
increases. 

In order to derive the form for the randomly-scattered superposed 
rays, mentioned in I) above, the difhise reflectors are assumed to be 
homogeneously distributed since they can be viewed as scattering the 

20 rays from randomly placed locations. Moreover, at any given path length, 
all ray amplitudes are approximately equal, and so the random scattering 
is only manifested In the phases of individual received rays. 

As time increases, the likelihood of any two received rays being in 
phase with one another is reduced, as rays with longer reflected paths are 

25 scattered further. Hence, the phases of individual rays of the same path 
length become increasingly divergent, and consequentiy, the amplitudes of 
the individual rays are more likely to cancel one another out upon 
superposition. 

The exponential decay of the random-scattering effect is governed 
30 by the equation: 
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m^Ae-"^' , (6) 
Where / is the propagation delay of the rays. The constant b 
determines how quickly the function tends to zero. Since more re-radiation 
occurs at shorter wavelengths, the wavelength, X , has the property: 

S feoc A. 

Replacing b with the physical parameter A yields the following form 
for the random phase approximation of the received signal: 

/(0=^e-*''^ (7) 
where A is the amplitude of the transmitted signal. 
10 In order to determine the number of dtffusers with respect to time, 

mentioned in ii) above, it is assumed that each diffuser is a point 
homogeneously distributed in space. For any fixed path length for a signal 
transmitted from a transmitter, ^ , to a receiver, B , the reflectors must lie 
on the surface of an ellipsoid. Therefore, the number of reflectors N is 
15 proportional to the surface area of an ellipsoid with foci separated by the 
distance AB. 

Figure 4 shows the surfaces that diffusers lie on for a fixed path 
length r , where rays travel from point A to point B. As r increases, the 
20 path from ^ to ^ becomes more haphazard due to the increased number 
of diffusers encountered. 

An analytical expression for the surface area of an ellipsoid exists, 
but contains a complete elliptic integral of the second kind, which has to be 
evaluated numericaiiy. However, an approximate analytical form for the 
25 surtece area can be derived, having a more intuitive form. The 
approximation for the surface area of the ellipsoid must hold in the limits 
t^cAB and r-><». 

In the case of ttie surfece area for the shortest reflection delay, 
which can equivalentiy be viewed as the surface area for AB-^oo^ the 
30 reflectors lie on the surface of a cylinder of area: 

S = 27AAB , (8) 
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Where b is the minor axis of the ellipsoid. 

Given that r is the reflected path length of a ray. then: 

^ 2 • (9) 

6 Therefore, for / ^ -> o , the number of reflectors is given by: 

In the case of the surface area for the longest reflection delay, 
which can equivalently be viewed as the surface area for ^-»o. the 
10 reflectors lie on the surface of a sphere, which has the simple form: 

Where the path length r is twice the radius of the sphere. 
Hence, the approximate form for the number of diffusers as a 
function of time Is given by: 

JV(0=^P^^p^+/^j . (12) 

where c is the speed of light. / is the duration for the received 
reflected rays to travel from A to B and is the time for a ray to travel 
directly from A to B. 
20 So, over time the number of diffusers increases virtien the surface 

they lie on evolves from a cylinder into a sphere. Thus, the number of 
received rays in time, and therefore the received signal, also increase by 
the factor N(t) due to the increased number of diffusers encountered by 
the rays. 

25 However, as time increases the amplitudes decay in free space, as 

mentioned In iii) above, by the amount: 

• (13) 
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The effect of N(t) is reduced by free-space propagation decay to give the 
amount 



(14) 



The effect of the three factors. 0 to ill), of random-phase 
supeiposition, the number of dlffusers, and free^pace propagation decay, 
thereby combine to yield the following approximation for the diffuse 
components of the received signal: 



10 



25 



/(0 = 5c-« 
Where 



This can be simplified to: 



+ 1 



(16) 



=^ and or=^ . 



The diffuse approximation can be viewed as a background function 
15 for the received signal r(0 . virtiere specular reflections represented by the 
linear superposition signal basis can be superposed, resulting in the 
transformation of the MEDLL equation into the form: 



m=0 



(17) 



As the diffuse scattering process is random, involves a large number of 
20 rays and exhibits an exponential decay, tiie above approximation provides 
a good characterisation of the diffuse process. 

The maximum likelihood estimator is given hereinbefore In Equation 

(1). 

The present invention hereby provides an Improved estimator, witii 
ttie capability to clean the correlation function of diffuse reflections, given 
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by the following equatlon: 




(18) 



5 Where s(0 is the replica signal, and B and a are free parameters. The 

other parameters all appear In the maximum likelihood estimator described 
hereinbefore. 



deriving from the approximation for the diffuse components of the received 
signal given in Equation (16), represents an improvement to the maximum 
likelihood (MEDLL) estimator. 

Figure 3 illustrates an example of the method by which the present 
15 invention may be implemented. 

It will be apparent to a person skilled in the art that this method 
could be embodied in either software or hardware, for example as part of 
the signal-processing unit of a radio receiver. 



The term 



10 




20 



